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Abstract 

The decay width, forward-backward asymmetry and lepton longitudinal and transver- 
sal polarization for the exclusive B — ► K*i + l~dec&y in a two Higgs doublet model are 
computed. It is shown that all these quantities are very effective tools for establishing 
new physics. 
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1 Introduction 



Experimental discovery of the inclusive and exclusive B — > X s 7 and B — > K*^ |T[] decays 
stimulated the study of rare B meson decays in a new manner. These decays take place via 
flavor- changing neutral current (FCNC) b — > s transitions which are absent in the Standard 
Model (SM) at tree level and appear only at the loop level. Therefore the study of these 
decays can provide sensitive tests for investigating the structure of the SM at the loop level, 
searching for new physics beyond the SM . 

Currently, the main interest is focused on the rare B meson decays for which the SM 
predicts large branching ratios and that can be potentially measurable in the near future. 
The rare B — > K*£ + £~(£ = e, fi, r) decays are such decays. For these decays the experi- 
mental situation is very promising || with e + e~ and hadron colliders focusing only on the 
observation of exclusive modes with I = e, fi and r as the final states. At the quark level, 
the process B — > K*£ + £~ is described by the b — > sl + l~ transition. In literature [4-17] 
this transition has been investigated extensively in both SM and two Higgs doublet model 
(2HDM). It is well known that in the 2HDM, the up type quarks acquire their masses from 
Yukawa couplings to the Higgs doublet H 2 (with the vacuum expectation value t> 2 ) and 
down type quarks and leptons acquire their masses from Yukawa couplings to the other 
Higgs doublet Hi (with the expectation value V\). In 2HDM there exist five physical Higgs 
fields: neutral scalar if , h°, neutral pseudoscalar A and charged Higgs bosons H ± . Such 



a model occurs as a natural feature of the supersymmetric models [IS]. In these models 



the interaction vertex of the Higgs boson and fermions depends on the ratio tan (3 = jjj 
which is a free parameter in the model. The constraints on tan ft are usually obtained from 
B — B, K — K mixing, b — > S7 decay width, semileptonic decay b — > ctv t and is given by 



la 20 



0. 7 < tan/J < .6(^_) (1) 

(the lower bound m H + > 200 GeV is obtained in |20|). 

In 2HDM the charged and neutral Higgs interactions with fermions induce new contri- 
butions to the FCNC processes. For b — > s£ + £~ {£ = e, jj) decay the contributions from 
neutral Higgs boson exchange are usually neglected, due to the fact that the interaction of 
the neutral Higgs and leptons is proportional to the lepton mass. But for the b — > st + t~ 
decay the mass of the r lepton is not too small compared to the mass of the 6-quark, 
therefore the neutral Higgs boson exchange diagrams can give considerable contributions 
to the process. Recently it has been emphasized by Hewett [^T[] that r lepton longitudinal 



polarization is an important observable and may be accessible in the B — > X s t + t~ decay. 
In |22| it is shown that in addition to longitudinal component Pl of polarization of the r 
lepton two other orthogonal components Pt, which lies in the decay plane, and Pn which is 
perpendicular to the decay plane, are also very significant for the r + r~ channel, since they 
are proportional to the mass of the r lepton. The polarization components Pl, Pt and P^ 
involve different combinations of the Wilson coefficients C7, Cg^ and C10 (see below) and 
hence contain independent information. Therefore, polarization effects can play a central 
role for the investigation of the structure of the SM and for establishing new physics beyond 
it. 



1 



In calculating the branching ratios and other observables at hadronic level, e.g., for 
B — > K*£ + £~decay, we face the problem of computing the matrix element of the effective 
Hamiltonian responsible for this decay, between B and K* states. This problem is related 
to the non-perturbative sector of QCD and it can be solved only by means of a non- 
perturbative approach. 

These matrix elements have been investigated in the framework of different approaches 
such as chiral theory |23| , three point QCD sum rules method 0] , relativistic quark model 
by the light-front formalism [p5|j , effective heavy quark theory |26[ and light cone QCD 
sum rules [f27| . The aim of the present work is to calculate these matrix elements using 
the light cone QCD sum rules formalism in the framework of the 2HDM, taking into ac- 
count the newly appearing operators due to the neutral Higgs exchange diagrams, and 
to study the forward-backward asymmetry and final lepton polarization for the exclusive 
B — > K*£ + £~decay. Taking into account the additional neutral Higgs boson exchange 
diagrams, the effective Hamiltonian is calculated in [28| as 
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v*v; a \ y Ec i (ii)o i (iM) + Y,CQ i (j*)Qi(t*) 



(2) 



where the first set of operators in the curly brackets describe the effective Hamiltonian 
responsible for the b — > sl + l~ decay in the SM. The value of the corresponding Wilson 
coefficients includes the contribution of the diagrams with H ± running in the loop (see f7j), 
i.e. 
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The explicit forms of C$ M (M W ), C$ M (M W ), and Cf^(M w ) and the corresponding oper- 
ators can be found in [^|. 

The second set of operators in the brackets, whose explicit forms are presented in p8[ , 
come from the exchange of the neutral Higgs bosons. The corresponding Wilson coefficients 
are: 
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The QCD correction to the Wilson coefficients Ci(mw) and Cq^ttiw) can be calculated 
using the renormalization group equations. In [^H] it was shown that the operators Og and 
Oio do not mix with Q i (i — 1, . . . , 10), so that the Wilson coefficients C 9 and C w remain 
unchanged and their values are the same as in the SM. Their explicit forms can be found 
in 1 28], where it is also shown that O7 can mix with Qj. But additional terms due to this 
mixing can safely be neglected since the corrections to the SM value of C7 arising from 
these terms are less than 5% when tan j3 < 50. 

Moreover the operators Oj (i = 1, . . . , 10) and Qi (i = 3, . . . , 10) do not mix with Qi 
and Q2 and also there is no mixing between Qi and Q2- For this reason the evolutions 
of the coefficients Cq 1 and Cq 2 are controlled by the anomalous dimensions of Qi and Q2 
respectively: 

C Qi {m b ) = v' 7Q/Po C Qi (m w ) , i = 1, 2, 

where 7q = —4 is the anomalous dimension of the operator s^r |P9| . 

Neglecting the strange quark mass, the matrix element for b — > s£ + £~ decay is | 28fl : 
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where q 2 is the invariant dileptonic mass. Here the coefficient Cg^(fi,q 2 ) = Cg(//) + 
Y(fi, p 2 ), where the function Y contains the contributions from the one loop matrix ele- 



ment of the four-quark operators |7|, [JO], |31| . It is clear that the last two terms in eq. (8) can 
give considerable contribution only for the B —>■ K*T + T~decay mode, since Cq 1 and Cq 2 
are proportional to the lepton mass mg. In addition to the short distance contributions, it 
is possible to take into account the long distance effects associated with real cc in the inter- 
mediate states, i.e., the cascade process B — > K*J/ip(ip') — > K*£ + £~. These contributions 
are taken into account by introducing a Breit-Wigner form of the resonance propagator and 



this procedure leads to an additional contribution to Cp of the form [02 
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From eq.(8) it follows that, in order to calculate the branching ratio for the exclusive 
B — > K*T + T~decay, the matrix elements (K* |s7 M (l - 7 5 )6| B), (K* \sicr^ u q u (l +7 5 )&| B), 
and (K* + 7s)&| B) have to be calculated. These matrix elements can be written in 
terms of the form factors in the following way: 



(K*(p,e)\s llM (l- l5 )b\B(p + q)) 
2V(q 2 ^ 



m B + m K * 



- ie*(m B + m K *)Ax{q ) + i(2p + q) fl (e*q) 



2m 



K* 



A 3 (q 2 ) - A (q 2 



m B + m K * 



+ 



(9) 



(K*(p, e) Isia^il + 7B )6| B(p + q)) = 
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where e is the polarization vector of K* meson. 

To calculate the matrix element {K*\s(l + j 5 )b\B) we multiply both sides of eq.(9) by 
g M and use the equation of motion. Neglecting the mass of the strange quark, we get: 
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(11) 



Using the equation of motion, the formfactor A 3 (q 2 ) can be written as a linear combination 



of the formfactors Ai and A 2 (see [24 
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Making use of this relation we obtain: 



(K*(p, e) \s{l + 75 )6| B(p + q)) = ^ {-i(e* q) A (q 2 )} 



(13) 



From eqs.(8), (9), (10) and (13) we obtain for the matrix element of the B — > K*£ + £ decay: 
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The formfactors V, Ax, A 2 , A , T l5 T 2 and T 3 are calculated in the framework of light-cone 
QCD sum rules in p7| and their q 2 dependence, to a good accuracy, can be represented in 
the following pole form: 
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which we will use in further numerical analysis (for A% we use eq.(12)). 

Using eq.(14) and performing summation over final lepton polarization, we get for the 
double differential decay rate: 
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where z = cosd , 9 is the angle between the momentum of the I lepton and that of the B 
meson in the center of mass frame of the lepton pair, A = 1 + r 2 + s 2 — 2r — 2s — 2rs, 

r = m K* s — anc i v — i — me are the lepton velocity and mass respectively. 
m B m B V 1 
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If we put Fi = and F 2 = (which corresponds to the case where the contributions of 
the neutral Higgs boson exchange diagrams are neglected) we get the results of |27 ] . As we 
noted earlier, the forward-backward asymmetry and final lepton polarization asymmetry 
are very useful tools for extracting more precise information on the Wilson coefficients 
C7, Cg and Cio and also for searching new physics. Therefore in this work we shall study 
these quantities as well. The forward-backward asymmetry Afb is defined in the following 
way: 
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Let us now discuss the final lepton polarization. We define the following three orthogonal 
unit vectors: 
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where p\ and px* are the three momenta of the l~ and the K* meson, respectively, in the 
center of mass of the £ + l~ system. The differential decay rate for any given spin direction 
n of the l~ lepton, where n is a unit vector in the t~ lepton rest frame, can be written as 
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where the subscript "0" corresponds to the unpolarized case and it can be obtained from 
eq. (17) by integration over z. From definition of e*j it is obvious that Pt lies in the decay 
plane whose orientation is organized by the vectors p\ and px* , and Pn is perpendicular to 
this plane. 

The polarization components Pi (i = L, T, N) are defined as: 
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After standard calculations for P, (i = L, T) we get 
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The denominator A in eqs.(20) and (21) can be obtained from eq.(17) by integration over 
z of the terms within the curly bracket. We also calculate the Pjy component of the lepton 
polarization, but the numerical results show that its value is quite small and because of 
that we do not present its explicit form here. 

2 Numerical Analysis 

The values of the main input parameters, which appear in the expression for the decay 
width are: m b = 4.8 GeV, m c = 1.35 GeV, m T = 1.78 GeV, m M = 0.105 GeV, A QCD = 
225 MeV, rris = 5.28 GeV, and = 0.892 GeV. We use the pole form of the formfactors 
given in eq.(16). For B meson lifetime we take t(B ( i) = 1.56 x 10~ 12 s The values of 
the Wilson coefficients Cj M {m b ) and Cf M (mb) to the leading logarithmic approximation 
are ]34[ [35|| : 

C 7 = -0.315 , do = -4.642 . 

The expression Cg^ for the b —>■ s transition in the next to leading order approximation is 
given as (see for example PH) 

c 9 e// K) = 

Cl M {m b ) + Cf ~ (m b ) + 0.124iu(s) + g{m c , s) (3d + C 2 + 3C 3 + d + 3d + C 6 ) - 
- ±g(m q , s) (d + 3d) - ^g(m b , s) (4C 3 + 4C 4 + 3d + C 6 ) + 

+ I (3d + d + 3C 5 + d) , (22) 
with 

d = -0.249 , d = 1-108, d = 1-112 x 10" 2 , d = -2.569 x 10~ 2 , 
d = 7.4 x 10" 3 , d = -3.144 x 10~ 2 , C^ M (m b ) = 4.227, 

Tfl f7 2 

where rh g = s = -±-7. 
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In the above expression, w(s) represents the one gluon correction to the matrix element 
Og and its explicit form can be found in [13[], while the function g(m q ,s) arises from the 



one loop contributions of the four quark operators 0% — Oq (see for example [[34], |55|), i.e. 

g 8 4 2 

g(m q ,s') = -ylnm q + — + -y q --(2 + y q )^ll-y q + 

+ - y q ) hn ±±^ZM - zttJ + 6(y q - l)arctan-j=L=^ , (23) 

where y Q = ^S 2 -, and s' = 

s' m b 

In all numerical calculations we use the following values for the masses of the Higgs 
particles: m H ± = 200 GeV, m h o = 80 GeV, m H o = 150 GeV, and m^o = 100 GeV. We also 

take sina = and for tan (3 we choose the following set of values: tan (3 = 1, tan (3 = 30 
and tan j3 = 50. 

In Fig.l we present the q 2 dependence of the differential branching ratio for — > 
K*t + t~ decay with and without the long distance effects. It follows from this figure that 
the differential branching ratio is sensitive to the value of tan (3 if the if}' mass region is 
excluded. For example at tan (3 = 50 the differential branching ratio is approximately 2 
times larger than the one at tan [3 = 1. 

In Fig. 2 we plot the dependence of the forward-backward asymmetry A FB on q 2 with 
and without the long-distance effects at different values of tan [3. From this figure we see 
that A FB is positive for all values of q 2 except in the if}' resonance region and it is sensitive 
to the value of tan (3. 

In Fig. 3 we depicted the q 2 dependence of the longitudinal polarization of the final 
lepton Pl with and without the long distance effects at different values of tan (3. As we see 
from this figure if we exclude the resonance mass region of if)' , Pl is negative for all values 
of q 2 at tan (3 — 1 as well as at tan [3 = 30 and tan [3 = 50. 

In Fig. 4 we present the q 2 dependence of the transversal polarization P F of the r lepton 
which lies in the decay plane, without long distance effects at tan (3 = 1, tan (3 = 30 and 
at tan (3 = 50. From this figure it follows that at tan (3 = 1 P T is positive at all values 
of q 2 . For tan (3 = 30 and tan (3 = 50, Pp is positive near the threshold region but far 
from the threshold it becomes negative. Therefore the determination of the sign of Pp in 
future experiments is a very important issue and can provide a direct information for the 
establishment of new physics. 

In Fig. 5 we present the differential branching ratio versus q 2 for B — > K*fi + fi~. From 
this figure and Fig.l it follows that the differential branching ratio for the B — * K*fi + fi~ decay 
is approximately one order of magnitude greater then that of the B — > K*r + r^dec8ij. Note 
that for tan (3 = 30 and tan (3 = 50, the results for the differential branching ratio are 
practically the same. 

In Fig. 6, we depict the dependence of A FB (B — > K*u + u~) on q 2 at tan [3 = 1 and 
tan (3 = 50. From a comparison of this figure with Fig.2 we observe that the behaviour of 
Afb for the B — > K*{x + u~dec&y is totaly different from that for the B — > K*r + T~dec&y. In 
the B — > K*T + T~case A FB is positive for all values of q 2 (without the long distance effects), 
while it changes its sign here in the B — > K* u + u~case. This indicates that for the B — > 
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K*r + T~case, the neutral Higgs boson exchange diagrams give considerable contribution. It 
follows from these two figures then that the determination of the sign of Afb in different 
kinematical regions is also a very important tool for establishing new physics. 

The behaviour of the longitudinal polarization with changing q 2 is presented in Fig. 7. 
In B — > K* fi + decay, Pl is firstly positive and then it changes sign (without the long 
distance effects), which is absolutely different from that of the B — » K*t + t~ case, for which 
Pl is negative for all q 2 . 

Since, as we have already noted, A FB , P L , and Pt contain independent information, 
their investigation in experiments will be a very efficient tool in establishing new physics. 

At the end of this section we present the values of the branching ratios for the Bd — > 
K*t + t~ decay. After integrating over q 2 we get for the branching ratios for the B& — > 
K*t + t~ decay, without the long distance contributions, 



B{B a 



K* 



T T 



0.86 x 10~ 7 (tan (3=1) 
< 0.90 x 10~ 7 (tan (3 = 30) 
1.11 xlO" 7 (tan ,3 = 50) 



(24) 



The ratio of the exclusive and inclusive channels is defined as 

B(B d -> K*t+t-) 



R 



B(b 



ST T 



-7 



In the SM this ratio is given as R = 0.270 ±0.007 when B(b -> sr+r") = (2.6 ±0.5) x 10 
35fl . In our case we have 

0.33 (tan (3=1) 
R = < 0.35 (tan (3 = 30) (25) 

0.43 (tan (3 = 50) , 

where we use the SM value for the inclusive B(b —>■ st + t~). 

In conclusion, we calculate the rare B K*£ + £~ decay in 2HDM. It is observed that 
Afb, Pl and the transversal polarization Pt of the charged lepton are very sensitive to the 
value of tan (3. Therefore, in search of new physics, their experimental investigation can 
serve as a crucial test. 
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Figure Captions 



In Fig.l, Fig.2 and Fig.3, lines 1, 3 and 4 correspond to the short distance contributions 
for the values of tan (3 = 50, tan (3 = 30 and tanf3 = 1, respectively; while lines 2, 5 
and 6 correspond to the sum of the short and long distance contributions for the values of 
tan f3 = 50, tan (5 = 30 and tan (3 = 1, respectively. 

1. Invariant mass squared (q 2 ) distribution of the r lepton pair for the decay B — * K*t + t~. 

2. The dependence of the forward-backward asymmetry A FB on q 2 for the decay B — > 

K*T + T~. 

3. The dependence of the longitudinal polarization P L on q 2 for the B — > K*t + t~ . 

4. The same as in Fig.3 but for the transversal polarization P T . 

5. The same as in Fig.l, but for the B — > K*\i + \T decay. 

In Fig. 6 and Fig. 7, lines 1 and 3 correspond to the short distance contributions for the 
values of tan (3 = 50 and tan (3 = 1, respectively; while lines 2 and 4 correspond to the sum 
of the short and long distance contributions for the values of tan (3 = 50 and tan (3 = 1, 
respectively. 

6. The same as in Fig.2, but for the B — > K*fi + fi~ decay. 

7. The same as in Fig.3, but for the B — > K*fj, + fi~ decay. 
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